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Abstract 

We investigate boundary multifractality of critical wave functions at the Anderson metal-insulator transition in two- 
dimensional disordered non-interacting electron systems with spin-orbit scattering. We show numerically that multifractal 
exponents at a corner with an opening angle 9 = 37r/2 are directly related to those near a straight boundary in the way 
dictated by conformal symmetry. This result extends our previous numerical results on corner multifractality obtained for 
9 < TT to > TT, and gives further supporting evidence for conformal invariance at criticality. We also propose a refinement 
of the validity of the symmetry relation of A. D. Mirlin et ai, Phys. Rev. Lett. 97 (2006) 046803, for corners. 
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Anderson metal-insulator transitions are continu- 
ous phase transitions driven by disorder. Examples of 
localization-delocalization (Anderson) transitions oc- 
curring in two-dimensions (2D) include non-interacting 
electronic systems with spin-orbit scattering ('sym- 
plectic symmetry class'), with sublattice symmetry, or 
in strong magnetic fields (quantum Hall effect). 

Recently, we have reported numerical evidence for 
the presence of conformal invariance at the 2D Ander- 
son transition in the symplectic symmetry class [1] . To 
that end, we have considered multifractal properties of 
critical wave functions near boundaries of disordered 
samples of finite size, and verified numerically that 
the multifractal exponents of critical wave functions 
at corners with opening angle (corner multifractal- 
ity) are related, through simple relations derived from 
conformal invariance, to the exponents computed near 



straight edges (surface multifractality). In Ref. [1] we 
have discussed corner multifractality at wedges with 
angles < n only, both acute {9 — 7r/4) and obtuse 
{9 = 37r/4). In this paper we extend this analysis to a 
corner with = 3tt/2 to show that the same equation 
relating surface and corner multifractality holds for the 
corner with a reflex angle {9 > n). 

Following Refs. [1,2], we define bulk, surface, and 
corner multifractality from the scaling of moments of 
wave functions tp{r) in bulk (b), surface (s), and corner 
(9) regions, 

L-'-W^^L-^^, (x = 0,s,b), (1) 

where dx is the spatial dimension of each region (db ~ 
2, ds = 1, and de = 0). The overbar represents the 
ensemble (disorder) average and the simultaneous spa- 
tial average over a region x surrounding the point r. 
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The exponents r^, r^, and are the bulk, surface, 
and corner muhifractal exponents, respectively. From 
the multifractal exponents we extract non-vanishing 
anomalous dimensions A^, 



2q + dx. 



(2) 



The multihactal singularity spectra /"(a) are obtained 
from Tg by Legendre transformation. 



/X / X\ X X 

[a ) q-Tg 



(3) 



As explained in Ref. [1] , under the assumption that 
the g-th moment |i/)(r)P'' is represented by a primary 
operator in an underlying conformal field theory [3], 
one can derive, using the conformal mapping w = z^^^ , 
the relation between the surface and corner multifrac- 
tal spectra /"(a^), 



7r , ^ 

■2 = ^K 



-2), f{al)^^[r{al)~l]. (4) 



The validity of these relations provides direct evidence 
for conformal invariance at a 2D Anderson transition 
and for the primary nature of the operator. 

In Ref. [1] we have shown that the probability dis- 
tribution of In |?/j(r)|^ becomes broader, as the opening 
angle 6 is reduced. This implies that the distribution 
is narrower at a corner with larger 9. We may thus 
expect that multifractal exponents can be more accu- 
rately calculated for corners with reflex angles than for 
corners with angles 9 < tt. We can then estimate the 
surface f^{aq) by taking the numerical data for and 
/^(Oq) obtained for 9 > n as input into Eq. (4). More- 
over, we can relate multifractal spectra of corners with 
different angles 9 and 9' {9 < 9'), yielding 



-(q, - 2), / (a,) 



:f\4)- (5) 



As we pointed out in Ref. [1], Eqs. (4) and (5) are 
valid only if all occurring ^ 0, because a, is non- 
negative for normalized wave functions. Thus, when 
the prefactor 9' 19 is larger than one (and hence ^ 
< ), the first of Eq. (5) cannot be used for q > 
g^, where q^ is a solution to = in Eq. (5). (We do 
not know if is finite for ^ ^ tt.) Taking this physical 
constraint into account, we find the following relation 
between anomalous dimensions for corners (9 < 9'), 



.6' 



—A" 



2{q 



1 < qe, 



q > 



(6) 



If we set 9 = n in Eq. (6) , we obtain a relation between 
anomalous dimensions at a surface {9 — tt) and a cor- 
ner with a refiex angle {9' > tt). In Ref. [1] we also dis- 
cussed the symmetry relation of Ref. [4], A, — AJ_q, 
and its application to corners x = 9. Here we propose, 
as a refinement of that discussion, that this symmetry 
relation (i) is valid for corners of any angle 9 including 
9 = n (straight boundaries) , but only in the range of q 
satisfying l — qg^q^qg, corresponding precisely [1,4] 
to the range ^ Qg ^ 4, and (ii) makes no statements 
about Aq for values of q outside of this range. [The de- 
pendence on q of Aq is linear for q > qg (corresponding 
to the termination of the multifractal spectrum [5]), 
whilst it may, in general, continue to be non-linear for 
q < I — qg, even [1,4] when a* > 4.] 

In this work, we numerically verify these relations by 
computing corner multifractal spectra at 9 = 37r/2 for 
the L-shape samples shown in the inset of Fig. 1. We 
take a tight-binding model with both random on-site 
potential and random SU(2) hopping [6], and numer- 
ically obtain, with the forced oscillator method [7], a 
wave function tp having energy eigenvalue closest to a 
critical point Ec = 1.0 (in units of the mean hopping) 
for each random realization characterized by the on- 
site disorder strength Wc ~ 5.952. The system size L is 
varied through L = 24, 30, ■ ■ ■ , 120 and the number of 
disordered samples is 6 x 10"' for each L. We set w = 2 




Fig. 1. (color online) Multifractal spectra /(a) for corners 
with 9 = 37r/2 (circles) and 8 = it/2 (squares), and surface 
(triangles) regions. Error bars are plotted at integer values 
of q for the corner with 9 = 3n/2. The solid and short- 
-dashed curves represent the theoretical prediction from 
Eq. (4), where f^'^^'^(ag^^^) and fict^) are used as input, 
respectively. The dashed curve is calculated from Eq. (5) 
with f^'^^^{ctq^^^) used as input. Inset: L-forni geometry 
with 3L^/4 sites. The shaded part is the corner region with 
9 = 3n/2 of the size Sw^/A. 
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Fig. 2. (color online) The exponents Ag/[q{l — q)] for corner 
with 8 = 37r/2 (circles) and 6 = tt/2 (squares), and surface 
(triangles) regions. Solid and dashed curves represent the 
conformal relation (6). 

of the corner region shown in Fig. 1. Muhifractal spec- 
tra are computed in the same way as in Ref. [1] . 

In Fig. 1, we show multifractal spectra /(a) of cor- 
ners with 9 = 3tt/2, together with those of corners 
with 6 = n/2, and of the surface region [1]. The peak 
position of f''l'^{aq''^) is a^''^ = 2.265 ± 0.003, 
which is smaller than cvq''^ = 2.837 ± 0.003. Also, the 
width of /^''/^(a) is smaller than that of f{a\). This 
is consistent with Eq. (4) at ^ = 37r/2. Figure 1 clearly 
shows that /"^'^''^(a) computed directly for the cor- 
ner with Q — 37r/2 agrees well with the short-dashed 
curve obtained from Eq. (4) while using f^{a) as input, 
which verifies Eq. (4) derived from conformal invari- 
ance. We have also calculated the surface f^{aq) and 
corner /^^'^{ciq^'^) from Eqs. (4) and (5), respectively, 
using f'^^^'^{ag^^'^) as input into these equations. This 
allows us to estimate f^ictq) near a'^ ~ 0, providing an 
estimate for = df^{a^ = 0) /da^. The theoretical pre- 
dictions (solid and dashed curves) are in good agree- 
ment with the numerical data (triangles and squares) 
for f^iccq) and f^^'^(aq^'^), respectively. 

Figure 2 shows the anomalous dimensions for 
corners with 6 — 3tt/2 and 9 = n/2, and the surface 
region, which are numerically calculated from the scal- 
ing |V(r)|29/(|^(r)|2)9 L"-^?. The solid and dashed 
curves represent the theoretical prediction, Eq. (6), 
from the conformal mapping using Aq^^^ as inputs. 
The data points for (triangles) agree with the solid 
curve for \q\ < 1.5, while those for A,''^ (squares) are 
close to the dashed curve only near g ~ 0. The data 
points for Aq^^^ satisfy, within error bars, the symme- 
try relation [4] Aq^^^ = A^^^q in the vicinity of q = 



1/2, indicating good numerical accuracy. This opens 
the possibility that one can use corner multifractality 
for ^ > TT to obtain, with the help of Eqs. (4) and (6), 
more accurate estimates for multifractal properties at 
a straight surface and corners with 9 < n. 

We briefly comment on multifractality of a whole 
sample with boundaries. We have found in Refs. [1,2] 
that corner multifractality may dominate multifractal- 
ity of a whole system, even in the thermodynamic limit, 
for large values of |q| if < Tq,Tq. Here we point 
out that this cannot happen with corners of reflex an- 
gles {6 > tt). The proof goes as follows. We first note 
that, from Eqs. (2) and (6), the difference of corner and 
surface multifractal exponents is given by — r| = 
1 -I- A^(7r/6' - 1) as long as > (note that > 
for 9 > n). Thus, when n < 9 < 2n, the inequality 
Tq > Tq holds If A^ ^ 2. Secoudly, since Tq is a con- 
vex function of q with the constraints tq = — dx and 
rf = 2 — dx (recalling [1] fi = 0, and thus A^ = 0), 
we find A^ s; for \q - 1/2| > 1/2 and < A^ < 2 
for < g < 1. We thus conclude that > r| when 
aq > 0. Finally, when Tq — Tq is positive for > the 
difference remains positive even if < oo and in the 
regime q ^ q^ where — 0, because r| is then con- 
stant for q ^ qs (and dTq/dq = ciq ^ aq). Hence, con- 
tributions from corner multifractality at 6^ > tt cannot 
be larger than contributions from surface multifractal- 
ity. The numerical results shown in Fig. 1 are consis- 
tent with the above general argument. 

In summary, we have investigated corner multifrac- 
tality for the reflex angle 6^ = 371/2 and confirmed the 
validity of the conformal symmetry relations. This re- 
sult provides stronger evidence for the presence of con- 
formal symmetry at the 2D Anderson metal-insulator 
transition with spin-orbit scattering. 
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